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Abstract 

A two dimensional spring-block type model is used to model capillarity driven self-organization of nanobristles. The 
model reveals the role of capillarity and van der Waals forces in the pattern formation mechanism. By taking into ac- 
count the relevant interactions several type of experimentally observed patterns are qualitatively well reproduced. The 
model offers the possibility to generate on computer novel nanobristle based structures, offering hints for designing 
further experiments. 
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Reproducible nanoscale patterns and structures are of 
wide interest nowadays for engineering components in 
modern small-scale electronic, optical and magnetic de- 
vices HI]. The so-called bottom up approach for the 
fabrication of these nanostructures uses nanoparticles as 
elementary building blocks. Under some specific con- 
ditions the nanoparticles self-organize into the desired 
structures |2|]. A well-known and widely explored pos- 
sibility to induce this self-organization is to use the cap- 
illarity forces which appear during the drying of a liquid 
suspension of nanoparticles Ollfjl]. For instance, regular 
and irregular two-dimensional polystyrene nanosphere 
arrays on silica substrates are generated by such meth- 
ods J5[]. These patterns are used then as a convenient 
mask in the NanoSphere Litography (NSL) method. 

Carbon nanotubes (CNT) are attractive materials for 
nanotechnology because of their interesting physico- 
chemical properties and molecular symmetries. In or- 
der to make them appropriate for certain applications, 
proper initial CNT configurations have to be built, and 
specific conditions have to be found which enable their 
controlled self-organization 0, [H . This is a very am- 
bitious and challenging task, which can be made easier 
by elaborating working computer models for the self- 
organization of CNTs on substrates. Therefore, not only 
experimental, but also computational studies can ad- 
vance the field of nanoengineering. 

In the work of Chakrapani et al. |@1 an experi- 
mental procedure is presented in which capillary self- 
organization of CNTs leads to puzzling cellular pat- 




Figure 1 : Scanning electron microscope images of the structures ob- 
tained in drying nanobristles. 



terns. After the collapse and reassembly of the highly 
ordered, anisotropic, elastic CNTs, shrinkage and crack 
formation occurs. The resulting cellular foams are visu- 
ally striking, stable patterns. 

The experimental procedure jgj] may be shortly sum- 
marized as follows. Multi-walled nanotube arrays are 
grown on rigid silica surface by chemical vapor depo- 
sition (CVD) based on the decomposition of ferrocene 
and xylene. The resulting nanotubes have a wall thick- 
ness of ca. lOnm and a diameter of ca. 30nm. The 
average distance between two nanotubes is ca. 50nm. 
The obtained nanotube bristle is oxidized in an oxygen 
plasma at room temperature and 133 Pa pressure for 10 
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minutes and immersed in a wetting fluid. After the liq- 
uid evaporates, characteristic cellular type patterns are 
formed, i.e. the ends of nanotubes self-organize in com- 
pact walls. 

Figure Q] shows scanning electron microscope im- 
ages of some typical structures. From the figures we 
learn that a wide variety of structures are engineered 
in such manner. Both statistically symmetric polygo- 
nal cells and rather elongated ones can be obtained by 
chaning the experimental conditions. 

These micrometer scale structures have many advan- 
tageous features. They can be elastically deformed, 
transferred to other substrates or floated out to produce 
free-standing macroscopic fabrics. Thus, they might 
find potential applications as shock absorbent reinforce- 
ment in nanofiltration devices, elastic membranes and 
fabrics, and containers for storage or growth of biologi- 
cal cells. 

Despite of its applications and the existence of well 
elaborated production protocols, the exact mechanisms 
responsible for self-organization of CNTs into vertically 
aligned cellular structures is not clearly understood. 
Recently, it has been argued that although we lack 
some basic information regarding the self-organization 
of CNTs within a bristle, this process can be approxi- 
mated with the self-organization of arrays of CNT mi- 
cropillars of micron-scale diameters 1 1(111 each consist- 
ing of thousands of CNTs. This observation enables 
the construction of a computationally tractable model 
which operates instead of stand-alone CNTs with mi- 
cropillars. 

In the present work a simple mechanical spring-block 
type model defined at mesoscopic micropillar level 
is considered for understanding the capillarity driven 
self-organization of nanobristles (or so called "CNT 
forests"). The aim is to build a model that is able 
to qualitatively reproduce the variety of experimentally 
produced patterns. 

The model used here for describing the self- 
organization process of nanobristles is a mesoscopic 
one and is based on the mechanical spring-block stick- 
slip model family. This model family ap pea red in 
1967, when R. Burridge and L. Knopoff [11] con- 
structed a simple mechanical model for explaining the 
Guttenberg-Richter law for the distribution of earth- 
quakes after their magnitude. The basic elements of 
the model are blocks and springs that interconnect in 
a lattice-like topology. The blocks can slide with fric- 
tion on a planar surface. The original system introduced 
by Burridge and Knopoff (BK) is a one-dimensional 
model. It can be studied numerically and it exhibits self- 
organized criticality [ 12]. 







Figure 2: Main elements of the spring-block model. Panel (a) shows 
a schematic 3D representation of the nanobristle for the initial and 
a later state. Panel (b) illustrates the dynamics of the equivalent 2D 
model. 



The BK model gained new perspectives with the 
strong development of computers and computer sim- 
ulation methods. Variants of the BK model proved 
to be useful in describing complex phenomena where 
avalanche-like processes are present, pattern forma- 
tion phenomena and mesoscopic processes in solid-state 
physics or material sciences 

ill til- 

Recently, by using this model, we have succes- 
fully explained the patterns obtained in capillary self- 
organization of nanospheres JHQiJ. Motivated by this 
success, hereby we propose to map the capillarity driven 
self-organization of nanotube bristles to a spring-block 
system, and to understand the pattern selection process 
by means of computer simulations. 

First, let us consider the three-dimensional (3D) 
model, which is very similar to the real nanotube ar- 
rangement. As sketched in Figure [2f a), the micropil- 
lars composed by thousands of nanotubes having fixed 
bottom ends are modeled by flexible strands. Their 
capillary interactions are represented by non-classical 
springs that connect the neighboring pillars. The evapo- 
ration of the liquid is simulated by the stepwise increase 
of tension in the springs. This will result in the agglom- 
eration of micropillar ends creating the final structure in 
the studied system. 

As shown in Figure |2fb), this 3D model can be eas- 
ily mapped into a two-dimensional (2D) one by pro- 
jecting the micropillars' top ends on the surface. In 
the projection plane the micropillars bottom ends are 
represented as dots, and their positions are fixed on a 
predefined lattice. The movable top ends are modeled 
by the disk shaped blocks which can slide with fric- 
tion on the 2D simulation surface. For visual purposes 
only, each block is connected by an extensible string 
with its bottom end showing the micropillars' trunk. In 
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Figure 3: Forces acting between micropillars immersed into liquid 
(top panel), the length dependence of the special spring force used for 
modeling this interaction (solid line on bottom panel), and the hard- 
core repulsion force between blocks (dashed line on bottom panel) 



our simplest approach there is no restriction imposed 
to the length of these extensible strings which means 
that nanotubes with infinite length are used. This cor- 
responds to the real case when the nanotubes length is 
much grater than the linear size of the cells in the final 
patterns. The blocks (top of micropillars) are connected 
with their nearest neighbors through special springs that 
model the resulting forces acting between the micropil- 
lars. 

These special springs are one key ingredient of our 
computational model. They represent the resultant in- 
teraction force acting between two micropillars im- 
mersed in suspension. The tension force in the spring 
has a complex variation with the spring-length, it's form 
is sketched with a red line in the top panel of Figure [3] 
This force is the resultant of the capillary force, the elec- 
trostatic repulsion and the van der Waals attraction be- 
tween the micropillars. The form of the capillary force 
acting between two micro-scale rods wetted by a liq- 
uid was measured experimentally and deduced analyti- 
cally, too 117Lll8lll9ll . Electrostatic repulsion and van 
der Waals attraction are successfully described by the 



Derjaguin-Landau-Verwey-Overbeek theory [20]. The 
capillary force decays inversely proportionally with dis- 
tance F c cc 1 Id. The resultant of all these forces reduces 
to zero for long distances and it has a maximum F m at 
distance d m of micropillars. 

Our model approximates this resulting force. The 
special spring-force used in simulations is shown with 
solid line in the bottom panel of Figure [3] For small 
elongations d < d m this spring acts as a classical spring 
with 



F k {d) = k(d - 2R) , 



(1) 



where k is the spring constant and 2R is its equilibrium 
length. For elongations longer than d m the spring force 
decays as 



F k (d) = k'/d, 



(2) 



where the constant k' is selected in such way that the 
force-elongation curve is continuous at d,„. The maxi- 
mum of the spring force at d m is denoted by F m . 

Similarly with our previous models of drying gran- 
ular materials 1 1511 or self-organizing nanosphere sys- 
tems [16], the effects created by the evaporation of the 
liquid is introduced through these springs. As the liq- 
uid evaporates the meniscus accounting for the capillar- 
ity forces gets more accentuated. This is modeled by a 
step-by-step increasing of the spring constant k. It has 
to be noted that through the drying process the maxi- 
mum force in the spring F m remains unchanged since 
the electrostatic forces are not affected. Accordingly, 
with the increasing of k the d m value has to be propor- 
tionally lowered. 

The second key ingredient of our spring-block model 
(necessary to get realistic structures) is the capillary 
force resulting from the non-vertical orientation of mi- 
cropillars. Once the micropillar becomes inclined, the 
meniscus radius of the liquid surface in contact with the 
micropillar becomes grater on the top side than on the 
bottom side. This leads to an unstable state because a re- 
sulting net force F ~ l/cos8 pointing vertically down- 
wards will act at the crossing point between the liquid 
surface and the micropillar 12111 . tending to incline even 
more the tube. Here we denoted by 6 the angle between 
the vertical and the tangent to the micropillar at the liq- 
uid level. 

In our spring-block approach this force which is 
monotonically increasing with the inclination angle is 
approximated by a simple linear repulsion force 



Fa\?fy — k a x 



(3) 



3 



acting between the block and its fixed bottom end (x 
denotes their distance in the simulation plane and k a is 
the repulsion constant). 

There is an additional almost hard-core-type repul- 
sion Fj which forbids blocks to interpenetrate. This is 
taken into account by the repulsive part of a Lennard- 
Jones type force which acts only when the distance be- 
tween two blocks becomes smaller than 2R (dashed line 
on the bottom panel of Figure[3]l. 

Additionally to the presented forces, damping forces 
are considered to stabilize the dynamics. A friction (pin- 
ning) between blocks and surface is introduced. It can 
equilibrate a net force less than F p . Whenever the total 
force acting on a block exceeds F p , the block slips with 
an over-damped motion. 

The dynamics leading to pattern formation consists of 
the following relaxation steps. First, the system is ini- 
tialized. Blocks are placed on a triangular lattice with 
lattice constant a and their corresponding bottom ends 
are fixed at the same positions. Then, the blocks are 
slightly dislodged in a random direction with a small 
random shift not grater than the half of the empty space 
between blocks. Thereby, the initial imperfectness of 
the nanobristle is modeled. Next, the interconnecting 
spring-network is build. The spring constant k is se- 
lected in such way that in the initial system the spring 
forces are not exceeding the pinning force F p that acts 
on each block. An initially prestressed spring-block net- 
work is thus constructed. 

During each simulation step the spring constant is in- 
creased by a small amount 5k representing the increase 
of tension due to the evaporation of water and the sys- 
tem relaxes to an equilibrium configuration. In this con- 
figuration the total net force acting on each disk is lower 
in magnitude than the pinning threshold F p . 

Similarly with the modelling of drying nanosphere 
suspensions IU6I1 . the relaxation dynamics is realized 
through an over-damped molecular dynamics simula- 
tion using a fixed time-step. The equilibrium state in 
a viscous medium can be found by moving the blocks 
in each simulation step in the direction of the net force 
acting on them, and with a displacement which is pro- 
portional to the magnitude of the resultant force 

dr = -¥{r)dt, (4) 

v 

where v denotes a viscosity. 

Since the simulation steps will be chosen small 
enough, this simplified dynamics is able to replace the 
Newtonian solution without loss of significant informa- 
tion. We remind here that we are not interested in the 
dynamics, but in the final equilibrium configuration. A 



relaxation step is finished when no disk slipping event 
occurs for the given spring constant value. It usually 
takes a very long time to achieve a perfect relaxation, 
therefore we introduce a tolerance (10~ 6 -10~ 9 ), and as- 
sume that the relaxation is completed when the largest 
displacement per unit time is smaller than this value. 
After relaxation is done, we proceed to the next simu- 
lation step and increase all spring constants by 5k. This 
dynamics is repeated until a final, stable configuration 
of the blocks is reached. 

When implementing the above relaxation dynamics, 
several types of boundary conditions can be considered. 
However, as it was shown in lfl6ll the boundary con- 
ditions (free, fixed or periodic) will influence the final 
stable structure only in the vicinity of boundaries. In 
the bulk, the obtained structures are rather independent 
of this choice. Accordingly, in the present simulations 
fixed boundary conditions are used and snapshots from 
the bulk are taken for later investigations. 

The above presented model has many parameters. 
Here, their values for our present simulations are given. 
The disk shaped blocks are considered with radii R = 1 , 
and it defines the unit length in the system. After fix- 
ing the diameter of the circular simulation area D 
-100 - 600, the blocks are placed on the triangular lattice 
having a lattice constant a = 2.2 - 5.0. The density (or 
space filling) of micropillars is implicitly defined by this 
lattice constant. The block sliding dynamics is governed 
by the viscosity v = 250 used in the equation © and the 
pinning threshold F p = 0.001. The springs used in sim- 
ulations are characterized by two parameters, namely 
the initial spring constant k = 0.01 - 0.05 and the equi- 
librium distance of springs 2R. By these parameters the 
unit force in our model system is defined. In order to 
simulate a quasi-static drying process the spring con- 
stant increasing step has to be a small one 5k = 0.001. 
The central capillary repulsion constant k a = 0.0003 is 
set to be small enough that in the initial system, where 
the micropillars are almost vertical, not to affect the cell 
formation dynamics. Later, when the walls are forming 
(and the micropillars are no longer vertically aligned), 
this force becomes grater and it helps the wall forma- 
tion and stabilization process. The used Lennard-Jones 
force has its standard parameters set to <x = 1 .79 and 
e = 1.310 -7 expressed in simulation units. All results 
presented in the following are obtained with the above 
choosed parameter values, unless it is otherwise speci- 
fied in figure captions. 

The presented algorithm can be easily implemented 
and systems up to 300 000 nanotubes can be simulated 
in reasonable computational time. As observable from 
the time-sequence in Figure |4] the cellular patterns are 
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N = 15 



N = 16 



Figure 4: Time evolution of the simulation for parameters D = 400, 
a = 2.2 and k = 0.01. The simulation time step N is noted below the 
snapshots. 



formed after nucleation of voids in the spring-block net- 
work (time step N = 9). A preliminary void is enlarged 
by the tensioned springs (time step N = 13) until the 
top of micropillars arrange in a final and stable cellular 
structure (time step = 16). The obtained dynamics re- 
sembles the pattern formation mechanism known from 
experimental in-situ observations Jgt]. 

Furthermore, the effect of nanotube density on final 
patterns is computationally investigated. Simulations 
with the same parameter set are preformed for systems 
initialized with different lattice constants a. By this, the 
space filling of the micropillar system 



P 



2tt 



(5) 



is varied. Here, it has to be noted that the micropillar 
density and implicitly the micropillar lattice constant 
are linearly connected to the real nanotube density p„ 




Figure 5: Different final structures depending on micropillar density 
obtained by simulations with parameters D = 600 and k = 0.05. The 
lattice constant is a = 2.4 for the top left panel, a = 3.0 for the top 
right panel and a = 5.0 for the bottom left panel. The obtained struc- 
tures are magnified on the bottom right panel. 



and the nanotube lattice constant a„, respectively. If a 
micropillar is composed by N nanotubes, then by sim- 
ple geometrical calculations it can be shown that 



nR 2 



P=Pn 



2R 2 „^N 



and a — a„ 



(6) 



where R„ denotes the radius of a nanotube. 

In Figure[5]three different type of simulated structures 
are presented. For high space filling p = 0.688 corre- 
sponding to lattice constant a = 2.4 polygonal cellular 
structures are obtained similar to those on the left hand 
side of Figure [1] As one can observe on the magnified 
cell image, at the center there is a clean area formed by 
radially outgoing micropillars. For intermediate space 
filling p = 0.404 which corresponds to a lattice con- 
stant a = 3.0 the micropillars self-organize into elon- 
gated cellular structures. The obtained structures are in 
qualitative agreement with the experimental structures 
presented in the right hand side SEM images of Figure 
Q] When interpreting the image, one has to take into ac- 
count that the length scale of this snapshot is 1 .25 times 
smaller than the length scale of the previously discussed 
one. Accordingly, the typical size of the simulated elon- 
gated structures is 3 - 5 times grater than the typical 
size of the polygonal structures. For low space filling 
p = 0.145 corresponding to a lattice constant a = 5.0 
the micropillars form bundled clusters resembling the 
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Figure 6: Highly ordered experimentally designed patterns |23j] (top) 
in comparison with simulation results (bottom) for the parameter set 
D = 600, a = 2.6 and k = 0.01. 



recent experimental patterns of sequential assembly of 
self-similar nanotube clusters [22]. 

Finally our simulations explored also the possibility 
of designing highly ordered micropatterns in nanobris- 
tles. In this sense, an experimental procedure elabo- 
rated in the last years [23] has been reconstructed by our 
spring-block type simulations. The design procedure is 
based on the experimental observation that low-density 
regions or vacancies in the bristle play an important role 
in the pattern cell nucleation process. 

As shown in top panels of Figure |6l different kinds 
of highly ordered micropatterned structures have been 
created by etching regular vacancies with laser pulses. 
Capillary self-organization in such systems yields the 
structures presented on the top panels of Figure[6] 

Simulations on micropillar forests with vacancies of 
radius r = 5 created on rectangular and triangular lattice 
have been performed. The results are shown on bot- 
tom panels of Figure [6] In agreement with experimen- 
tal findings, our spring-block simulation results suggest 
that various micropatterns may be designed by proper 
preparation of the initial nanobristle taking into account 
that the wall of a polygon shaped cell forms approxi- 
mately at the vertical bisector of two adjacent vacancies. 

In conclusion, a simple mechanical, spring-block 
type model has been proposed here to model capillarity 
driven self-organization of nanobristles. The 3D prob- 
lem has been mapped to a 2D model that works at the 
mesoscopic micropillar level incorporating real interac- 
tions known from earlier experimental findings. Our 
computer simulations evidenced the role of capillary 
and electrostatic forces in forming the self-organized 
nanostructures. The dynamics leading to pattern forma- 



tion has been also revealed. By using the same model 
parameters with different nanotube densities three qual- 
itatively different types of patterns were reproduced. 
Moreover, the possibility of designing highly ordered 
nanostructures has also been explored. 
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